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^ , Abstract 

' The runaway breakdown - extensive atmospheric shower discharge (RB - EAS) 

excited in thunderstorm atmosphere by high energy cosmic ray particles (cp > 10^^ — 
10^^ eV) generate very powerful radio pulse. The RB - EAS theory is compared 
CSl ' with observations of radio pulses. An agreement between the theory and experiment 

^ ■ is established. The existence of nowaday satellite and ground based systems which 

ly-^ , obtain regularly a large amount of observational radio data could allow to use them 

I in combination with other methods for effective study of high energy cosmic ray 

2 ■ particles 

o 

(D : 1 Introduction 

Oh: 

1^ . Le Vine (1980) and Willet (1989) studying radio emission from thunderstorms have noted 

a distinct class of powerful radio pulses. The intensive investigations of this phenom- 
ena were performed during recent years. For the measurements the satellite FORTE 
^ I and specially constructed systems as LASA, EDOT and others were used [Smith et al., 

■ (1999), (2002); Thomas et al (2001); Light and Jacobson (2002); Jacobson (2003)]. These 

studies allowed to establish that the radio pulses having enormous power 100 - 300 GW 
are emitted in the wide frequency range by the intracloud discharges in the upper tropo- 
sphere. The pulses are short time (< lO/xs) and have a definite bipolar form. That is why 
they were called narrow bipolar pulses (NBP). The observations show that these strong 
radio pulses are accompanied by very weak optic emission only. A detailed analysis of the 
whole complex of observational data allowed Jacobson (2003) to state that NBP is a new 
type of thunderstorm discharge quite different from usual lightning. He speculated that 
it could have relevance to runaway breakdown effect. 

Runaway breakdown (RB) is a new physical concept of an avalanche type increase of 
a number of relativistic and thermal electrons in air proposed by Gurevich, Milikh and 
Roussel-Dupre (1992). The avalanche can grow in electric field E > E^. which is almost 
an order of magnitude less than the threshold of conventional breakdown. The electrons 
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with high energies e > {0.1 — l)MeV can become runaway and are accelerated under the 
action of electric field E > E^. Directly this process - acceleration and collisions with 
air molecules lead to the avalanche type growth of the number of runaway and thermal 
electrons (Gurevich and Zybin 2001). 

Runaway breakdown in air (RB) is stimulated by the presence of a high energy cosmic 
ray secondaries. Extensive atmospheric showers (EAS) are accompanied by an effective 
local growth of the number of cosmic ray secondaries and thus have a strong infiuence 
on the RB process (Gurevich et al 1999). The combined action of runaway breakdown 
and EAS lead to the development of RB - EAS discharge - new type of electric discharge 
where relativistic electrons play a decisive role. RB-EAS discharge is accompanied by 
strong exponential growth of the number of energetic and thermal electrons, positrons 
and gamma quants (Gurevich ct al 2004a). It can serve for the generation of a strong 
bipolar radio pulse (Gurevich et al, 2002). 

The goal of the present work is to extend a theory of RB - EAS discharge in air to a 
very high energy range of cosmic ray particles {sp > 10^^ — 10^^ eV) and to compare the 
theory with the results of NBP observations. 



Narrow bipolar pulses are isolated short time discharges generated in thunderclouds. They 
were discovered and first studied by Le Vine (1980) and Willett et al (1989). During last 
years very intensive and detailed studies of this special type of atmospheric discharge 
were performed by Smith et al (1999), Rison et al (1999), Thomas et al (2001), Light and 
Jacobson (2001), Smith et al (2003), Suszcynsky and Heavner (2003), Jacobson (2003) 
and others. It was established that NBP are compact and energetic intracloud discharge. 
They are observed in two forms: negative (NNBP) and positive (PNBP). For negative 
first electric field peak is negative and vice versa for positive. 
Main features 

1. NBP are high altitude discharges. Their main location (Smith et al 2003) 
z= 15-20 km for NNBP (sharp peak near 18 km) 

z= 7-15 km for PNBP (sharp peak near 13 km) 

2. Time characteristics of NBP (Smith et al 1999, 2002) 
mean rise time ~ 1 — 2//s 

full width at half maximum ~ (2 — 5)/is 
full rise + fall time ~ (5 — 10) fis 

3. NBP is observed as low frequency 0.2 - 0.5 MHz bipolar ground electromagnetic 
wave with large amplitude. The NBP effective field amplitude was measured simulta- 
neously at several stations situated at different distances R (Smith et al 2002, 2003). 
Characteristic wave field is 



4. Electric current pulse, generating NBP is unipolar and its maximum reaches values 



2 Narrow bipolar pulses (NBP) 




Jm ~ (30 - 100)M 
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5. Dipolc electric moment change in the cloud determined by NBP discharge is M 

0. 2 - 0.8 Cu km (Smith et al 1999). 

6. From the analysis of observational data it follows that the NBP current initiator 
moves with very high speed. The speed of initiator is determined as 7.3 x 10^ — 3.0 x 
10^°cm/s (Smith et al 1999) or ~ 10^°cm/s (Jacobson 2003). 

7. The NBP discharge rates is growing with the thunderstorm convection strength 
(Suszcynsky and Heavner 2003). 

HF radio emission 

1. NBP is always accompanied by intensive radio emission in a wide frequency range 
up to 500 MHz (Smith et al 1999). 

2. Detailed study of HF emission in the frequency range (26 - 48) MHz fulfilled at 
the FORTE satellite allowed to establish the following main features (Jacobson 2003): 

a. HF emission connected with NBP is always powerful its integrated ERP > 40 kW 
It is called by Jacobson strong intracloud (IC) pulse, 
b Strong IC pulses are incoherent in HF range. 

c. Strong IC pulses are accompanied by very low optic emission - at least two orders 
of magnitude less than in usual flashes. 

d. Strong IC pulses occur singly or initiate intracloud flashes. They never come within 
the interior of flashes. The optical emission does not occur for initiator of strong IC pulse. 

e. Strong IC pulse can occur without NBP, but not vice verse. 

f. The strong IC pulse in a given storm appear to have a truncated ERP distribution 
staying below a limiting ERP. For the most storms this is of the order of 1 MW in the 
FORTE HF band 26 - 48 MHz. The maximal value of a limiting ERP is 10 MW. 

3 RB - EAS discharge stimulated by high energy CR 
particles. 

1. RB - EAS discharge is determined by the exponential growth of the number of rela- 
tivistic electrons, positrons and gamma quants. The characteristic length is given in 
the Table for the main location heights of PNBP (z=13 km) and NNBP (z=18 km). In 

RB - EAS discharge the RB process is facilitated significantly in the energy range 3-30 
MeV due to effective generation of gamma quants and e+e" pairs (Dwyer 2003, Gurevich 
et al 2004a). 



Table: RB - EAS and NBP characteristics 



z{km) 


RB-EAS theory 


NBP observations 




la{m) Tg {lis) Trel {l^'S) 


Tg{lXs) TrelilJis) 


13 
18 


200 0.7 0.7 
400 1.4 3 


1-2 2-5 
1-2 2-5 



2. The production of a giant number of thermal electrons nth due to ionization of 
air molecules is going simultaneously with the generation of fast relativistic electrons 
nth ~ 10^^/) (Gurevich et al 2004). So the growth time of thermal and relativistic 
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electrons coincide. The attachment of thermal electrons to oxygen is due to three particle 
collisions. The characteristic attachment time (Phelps 1969): 




Tatt 

Here Nm - the number density of air molecules in cm Attachment determines the 
relaxation time of discharge. The growth Tg = la/c and relaxation time of RB - EAS 
discharge are presented in the Table. 

3. The electric current J generated in RB - EAS discharge is unipolar. For a given 
Em/ Ec its maximum Jm is proportional to the number of thermal electrons. The last is 
proportional to the number of fast electrons and thus proportional to the number of cosmic 
ray secondaries. That is why the maximal current Jm is proportional to the energy of 
cosmic ray particle Ep. For conditions of strong RB - EAS discharge, using the calculations 
performed in Gurevich et al 2004a, 2004b, we can estimate Jm as 



Jm-[-^]^^ (1) 



Here Sp is the energy of cosmic ray particle 

4. The radio pulse emitted by RB - EAS discharge is bipolar. Its amplitude depends 
on the distance to the source R and on discharge parameters. In average according to our 
simple model ((H) it is proportional roughly to Bp. 

5. The main energy dissipated in RB - EAS discharge goes to ionization of air molecules 
and to excitation of N2 vibration levels. Optic emission is very low - less that 1% of total 
energy (Roussel-Dupre and Gurevich 1996, Gurevich et al 2004). 



4 Discussion 

NBP and RB - EAS theory 

Let us compare NBP with the theory of RB - EAS discharge generated at high heights 
13 - 18 km by very energetic cosmic ray particles {Sp > 10^^ — 10^^ eV). Parameters of 
RB - EAS discharge at high heights are presented in the Table. One can see that both 
growth and relaxation time characteristics roughly agree with NBP observations. Electric 
current is unipolar. Its main peak reaches the values (P) which generally are in agreement 
with the average data obtained in NBP observations. One can see the following chain of 
measurements which show the rough proportionality between Jm and Ep-. 

Tien Shang experiment (Gurevich et al 2004b) 

ep~10^^eV, J^- 1-10 A 
Lightning first pulse measurements (Gurevich et al 2003) 

ep ~ 10^^ - lO^^eV, J^~0.1-lkA 

NBP 

- 10^^ - lO^^eV, J^~10-100kA 
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Thus RB - EAS discharge can serve as background for explanation of NBP phenomena. 

Note that the high values of maximal electric field E„J E^. ^ 1.2 — 1.4 were supposed in 
estimates Gurevich et al (2004a). The direct observations of electric fields at NBP heights 
gave lower values yet (Vonnegut et al 1989). On the other hand the NBP are usually seen 
in the active phase of storm near reflectivity core with 50 dBz radar reflectivity (Smith et 
al 1999). According to Smith et al 2003 the strong positively charged layer at the heights 
(15 - 16) km exist between generation maximum of NNBP and PNBP. It is impotant also 
that NBP formation is accompanied by powerful HF emission. 

HF emission model 

In Gurevich et al (2004a) RB - EAS model the air density distribution was sup- 
posed to be smooth function depending on one coordinate z. In reality in the high 
reflectivity core where NBP is generated a large number of hydrometeors (liquid or frozen 
water drops) exist. That can generate significant fiuctuations in RB process. Really ac- 
cording to MacGorman and Rust (1998) the characteristic hydrometeor dimensions in a 
high reflectivity core tq ~ 0.3 mm and density number n ~ 10~^ cm Due to this the 
mean free path of fast electrons between collisions with hydrometeors is l^h ~ (vrrgn)"^ ~ 
0.3 km. The effective radiation length of hydrometeors ^ O.Odg/cm'^ is equivalent to 
the length = th/ pair ~ 1.5 — 3 m in the air at the heights 13 - 18 km. It means that 
the passage of hydrometeor by runaway electron is equivalent to the passage of additional 
air length 1^. Thus the existence of randomly distributed hydrometeors lead to random 
fluctuations of and Em/Ec of the order of Ih/hh ~ l9o- Taking into account expo- 
nential dependence on Em/ E,. of fast and thermal electrons generated in RB process, one 
can expect that hydrometeors could be the reason of strong enough random fluctuations 
of electrons and electric current. 

Another effect which could be responsible for fluctuations is space inhomogeniety of 
cosmic ray secondaries which serve as a seeds of RB - EAS process. These fluctuations are 
especially significant at the high heights of the order 15-20 km where the shower has not 
reach the full stage of development yet. That is why the showers coming under inclination 
angle more than 60° are mostly effective The approximation of a quasi flat electromagnetic 
shower cascade moving together with cosmic ray particle which was used in Gurevich et 
al (2004a) could be considered as a good approximation for low energy Sp < 10^^ eV 
only. For the higher energies of cosmic ray particles situation is quite different. Note 
that laboratory experiments at accelerators in this energy region are absent. Both cosmic 
ray experiments and theory show very complicated character of shower structure at high 
energies Sp > lO^^eV (Murzin 1988). The observations using Pb and Pb-C cameras in 
Pamir and Chakaltaya experiments demonstrate multidimensional and inhomogeneous 
structures of electron - photon shower components (Baiburina et al 1984, Hasegawa and 
Tamoda 1996). Though photon super families of halo type (Genina et al 1981) and 
a number of new type processes "Centauro", "Chiron", penetrating chanels, anomalous 
transition curves (see Gladysz -Dziaoshu 2001, Capdevielle and Slavatinsky 1999 ) are not 
well understood, the growth of fluctuations in photon and secondary electrons distribution 
could exist. In spite of fluctuations the part of the energy of primary cosmic ray particle 
distributed in electron photon component of the shower is conserved at the level (0.2 — 
0.3)6p (Baburina et al 1984). It means that the full power of cosmic ray secondaries is 
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proportional to Sp what supports the basic relations (H)). 

Space and time fluctuations of Em/Ec and cosmic ray secondaries lead to the corre- 
sponding fluctuations in number of generated thermal electrons what results in a strong 
fluctuations of the current. In that a way we suppose the intensive IC pulse HF radiation 
is generated. 

At the same time these fluctuations do not change the essence of RB - EAS process. 
Its main features are reflected in IC HF radiation characteristics emphasized by Jacobson 
(2003). Radio emission is strong but the accompanying optic emission is very low, IC HF 
emission exists only as a single pulse (or as a first pulse in the beginnings of a flash ) and 
so on. 

Coherence effect and the giant power of low frequency radio emission pulse 
(NBP) 

The power dissipated in RB - EAS discharge is approximately proportional to the 
energy of cosmic ray particle Ep. The number of created newborn thermal electrons is also 
proportional to Ep and to the exponential factor F{Em/ Ec). Factor F in real conditions 
could reach values F ~ 10^ - 10^ (Gurevich et al 2003, 2004a). Thus the full energy Wd 
dissipated by thunderstorm electric field to the creation of thermal electrons and hence to 
the electric current of RB - EAS discharge is of the order of EpF. For Ep = 10^^ eV energy 
Wd ~ IJ, for Ep ~ 10^^ eV - Wd ~ 1 kJ. The same is fully correct for a non coherent HF 
radio emission. Its maximal frequency integrated power is 1-10 MW, maximal emitted 
energy (10 - 100) J. 

The picture changes dramatically when we take into account coherence of low frequency 
(LF) emission process. The pulse current region determined by the EAS scale is la ~ 
300 — 400 m. For LF radio emission (200 - 500) kHz the scale Lc is of the order or less 
than the length of radio wave A = (600 — 1500)m. It means that the LF radio emission 
is coherent process and due to this the power of radio emission in LF is growing with 
current J proportionally to J^. That is why the ground wave power in NBP could reach 
extreme values: 

2J2 

P = , P^ (100 - 300) GW 

3c 

It means that emitted by radio pulse energy can reach 0.2-1 MJ. We see that the work 
of thunderstorm electric field in NBP is really giant: it exceeds 10^ times the energy of 
cosmic ray particle triggering the process. 

We emphasize that the LF coherence effect is a strong argument supporting RB - EAS 
model of NBP generation. LF coherence means that though HF emission is incoherent 
but all currents which generate it according to the model has the same direction. In other 
words they have the same nature - thermal electrons moving by thundercloud electric 
field. These currents are only distributed in space - time inhomogeneously in groups 
and thin filaments - determined by inhomogeniety of local electric field and cosmic ray 
secondaries production. 

Conventional ohmic electric discharge in the night time ionosphere excited 
by NBP 

According to Smith et al (2003) in the night time conditions the NBP reflection region 
lies at the hieghts 2; 80 - 90 km. One can see, that NBP power is so high that electric 
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field of tlie pulse Enbp coud reach and overcome the value of conventional breakdown field 
Eth{z) at these heights 



where Nm{z) is the number density of air molecules at the height z. For the height 80 
km it follows that Eth = 24 V/m, for the height 90 km Eth = 3.2 V/m. Thus one can 
expect that short - time discharge (a few /xs) excited by NBP in lower ionosphere could 
be observed. 

Statistics of high energy cosmic rays and NBP events 

The flux of cosmic ray particles having the energy e > 5 x 10^^ eV is one per km^ per 
year (Murzin 1988). The annual thunderstorm duration complied from 450 air weather 
system in USA gives an average number 100 h/year (MacGorman et al 1984, Uman 
1987). It means that thunderstorms lasts about 1% of the whole time. The flux of high 
energy cosmic ray particles (e > 5 x lO^'^eV) in thunderstorm time is 0.01 /km"^ year. The 
FORTE satellite combined with ground based systems allows to detect simultaneously 
radio emission at (1 3)10^ km^. Thus the number of events of interaction of cosmic 
rays e > 5 x lO^^eV with thunderclouds is (1 3)10^ per year. In reality if the NBP 
can be generated by lower energies say 5 x 10^^ the statistics would be 100 times better 
~ 10^ events per year. Thus we see that there is no contradictions between the NBP 
observational data (10^ events during 4 years (Smith et al 2003)) and high energy cosmic 
ray flux. 

Initiator speed 

According to observations (Smith et all 1999, Jacobson 2003) the speed of NBP ini- 
tiator is close to lO^'' cm/s. The high energy cosmic ray particle (HECRP) moves with 
velocity of light. But the radio emission of RB -EAS discharge is generated by the current 
of the thermal electrons which flows in the direction of thunderstorm electric field E. Thus 
the speed of initiator v = ccos6, where 6 is the angle between E and HECRP direction 
of motion. That explains the difference in initiator speeds. We emphasize that values of 
observed velocity is not far from c. 

We note that for effective development of shower the cosmic ray particle has to pass a 
thick enough atmospheric layer. At the heights 13 - 18 km the density of air is low. That 
is why HECRP has to go under the large inclination angles 6 ^ 70° — 80° to obtain well 
developed shower. Thus the RB - EAS discharge become stronger for HECRP having 
large inclination angle. 

High energy cosmic ray particle detection (preliminary remarks) 

The high power of NBP signals allows to detect them at large distances - up to 1000 
km. That can be used for detection of cosmic ray particles with energies Ep ~ 10^*^ — 10^^ 
eV and higher (HECRP). Of course every concrete event depends not only on Sp but on 
a number of others factors like Em/Ec, reflectivity core size, thunderstorm activity. But 
one can expect that after averaging over a large number of events some useful information 
about HECRP could be obtained directly from NBP. For example, according to Smith 
et al (2003) the NBP database has now more 100 000 events. NBP is created by current 
having maximal value Jm which is proportional to (0). As the HECRP integral flux 
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F{e > Ep) is proportional to one can construct the averaged over all NBP events 
function 




Here $ is the full number of NBP with > Jq and Iq is arbitrary chosen current, say 5 
kA. A careful selection of PNBP and NNBP homogeneous conditions using HF emission 
data should be performed. At the middle values Jm < 100 kA the function Q{Jm) is 
expected to be approximately constant what will show that cosmic ray spectrum plays 
a dominant role. In that case the curve Q{Jm) and its behavior at the highest values 
Jm ^ 300 kA can give evidence of the existence of GZK cutoff or its absence. 

Here we supposed that up to highest values of flux dependence on Sp is decisive. 
The role of other thunderstorm factors should be carefully studied. 

5 Conclusion 

The analysis presented in the letter allows to formulate the following main statements: 

1. Narrow bipolar radio pulses (NBP) are generated by runaway breakdown - extensive 
atmospheric shower (RB -EAS) discharge initiated in atmosphere by very high energy 
cosmic ray particles. 

2. The detailed analysis of a large NBP data base could be used to obtain unic 
information about HECRP spectra including a highest energy range. There is a chance 
to establish the fundamental fact of existence (or non existence) of GZK cutoff. 

3. Of special interest is the detailed analysis of giant information gathered in HF 
radio spectrum observations. These new data analysis being combined with the theory 
and ground based measurements give a chance to reach some progress in understanding 
of basic processes in high energy particle physics {e > 10^^ eV). 
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